An exact renormalization group approach to frustrated magnets 



M. Tissier, B. Delamotte and D. Mouhanna 

Laboratoire de Physique Theorique et Hautes Energies. Universites Paris Vl-Pierre et Marie Curie - Paris VII-Denis 

Diderot, 2 Place Jussieu, 75252 Paris Cedex 05, France 

Frustrated magnets are a notorious example where usual perturbative methods fail. Having 
recourse to an exact renormalization group approach, one gets a coherent picture of the physics of 
Heisenberg frustrated magnets everywhere between d — 2 and d — 4: all known perturbative results 
are recovered in a single framework, their apparent conflict is explained while the description of the 
phase transition in d — 3 is found to be in good agreement with the experimental context. 



I. FRUSTRATED MAGNETS 



One of the great challenges of present-day condensed matter physics is to understand the effects of competing 
interactions - frustration - in three-dimensional spin systems. However, after twenty five years of investigations, the 
question of the nature of the universality class for the phase transition of the simplest frustrated modeL the classical 
antiferromagnetic Heisenberg model on a triangular lattice (AFHT model), is still strongly debateoo. The main 
feature of this system is that, due to frustration, its ground state is given by a configuration in which the spins have 
a non-collinear order: on an elementary triangular plaquette, they point at 120° one from another. A consequence of 
this non-collinear order is that, contrary to the case of non-frustrated spin systems, the microscopic order parameter 
is no longer a single vector but a matrix, built from two orthonormal vectors (^1,^2), playing a role similar to the 
staggered magnetization in non-frustrated antiferromagnetic systems. This implies a symmetry breaking scheme in 
which the group of spatial rotation is fully broken. In the case of the AFHT model it is given by 0(3)® 0(2) — > 0(2)cl. 
This situation raises the question of the nature of the phase transition for frustrated magnets in three dimensions. 
Experiments performed on materials supposed to belong to the AFHT universality class exhibit a critical behavior 
with exponents differing from those of the standard O(N) universality classu. However the experimental situation calls 
for several important remarks. First, the critical exponents differ significantly from one material to another and also 
with those obtained by Monte Carlo simulations performed either directly on the AFHT model or on related models 
(see ref.[4] for a review). Second, for a given material, the scaling relations are violated at least by two standard 
deviations. These results are manifestly not compatible with the existence of a true universality. 



II. PERTURBATIVE FIELD THEORETICAL CONTEXT 



From the theoretical point of view, this state of affairs raises the question of the ability of the field theoretical 
formalism to explain the experimental and numerical situation. As in the O(N) case, mainly two kinds of approaches 
have been used to tackle with this problem. We now display the results of these studies. 



A. The Non-Linear Sigma Model Approach 



The first approach consists in a low-temperature expansion performed around two dimensions on the suitable field 
theory, the 0(3)(£>0(2)/0(2) Non-Linear Sigma (NLct) modelB which is obtained by expressing the interaction between 
spins in terms of the Goldstone modes, the excitations associated to the transverse fluctuations of the magnetization. 
This is realized, as usual, by using the geometrical constraints of the ground state - here the orthonormality of <J)i 
and <f>2 - to get rid of the irrelevant - longitudinal - degrees of freedom. Formally, the resulting partition function 
writes: 



Z = I Dir exp 



-\J d d x gij (w) 3^3^ 



(1) 



where the tt 1, s represent the - three - Goldstone modes of the model, gij(ir) encoding their interaction. The expression 
(|l|) is then used to evaluate the renormalization of the interaction between the Goldstone modes in a double expansion 
in the temperature T and in e = d — 2. In the case of frustrated magnets, the result of this study is that there exists, 
above two dimensions, a stable fixed pointn. At this fixed point, the symmetry breaking scheme is enlarged from 
0(3) <g> 0(2)/ 0(2) to 0(3) <g> 0(3)/0(3). An important point is that the renormalization group (RG) properties of 
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a NLcr model depend only on the geometry of the order parameter space - the coset G/H - viewed as a manifold 
equipped with the metric gij(ir)u. In particular, the (3 functions for the coupling constants entering in the model are 
given in terms of geometrical quantities on G/H like Ricci - i£y - and Riemann - R? qr - tensors on the manifoldB: 

/%(flO = = e 9i] ~ ~ -^2 R i qT R ipqr ■ (2) 

This implies that the RG properties depend only on the local properties of G/H: the Lie algebras of the symmetry 
groups of the high and low temperature phases. As a consequence, they are insensitive to its global structure - 
topology - and to the representation used to realized the symmetry breaking scheme. It is known that the manifold 
0(3) ® 0(3)/0(3) has the same local structure as 0(4)/0(3) on which is built the NLcr model relevant to the four- 
component non-frustrated vector model. As seen from the NLcr model calculations, the two above theories are thus 
equivalent to all orders. Accordingly, the critical properties of the, the AFHT model between two and four dimensions 
should be governed by the standard O(4)/0(3) universality classB. This is however in strong disagreement with both 
experimental and numerical results as well as the weak-coupling analysis of the Landau- Ginzburg- Wilson (LGW) 
model of frustrated magnets around four dimensions. 

B. LGW model approach 

An alternative approach to that of the NLcr model is to consider a suitable LGW representation of the physics 
of frustrated magnets. This can be obtained, for instance, by relaxing, via a potential, the constraints defining the 
ground state: 




The P functions for the couplings u and v are-jdetermined in a double expansion in these couplings and in e = 4 — d. 
The result of this analysis, at dominant ordercTO, is that the RG equations obtained have no stable fixed point. 

C. A critical appraisal 

Frustrated magnets thus present a situation in which, contrary to the O(N) case, the different perturbative ap- 
proaches performed around the lower and upper critical dimensions conflict when naively extrapolated to three 
dimensions. Since there is no indication that any of these different perturbative results should fail in their respective 
domain of validity - i.e. for small e = d — 2 and small e = 4 — d - this situation raises two problems. First, a practical 
one, which is the determination of the critical behaviour in three dimensions. An high-order double expansion in 
the coupling constants and in e = 4 — d or a direct perturbative computation in three dimensions, which are both 
possible on the LGW model, can, a priori, deal with such a question. This depends, however, on the ability to extract 
reliable predictions out of the asymptotic series obtained, for example, by Borel resummation. But there is no general 
answer to this question. In fact, it now appears, from several recent studies, thatin the case of frustrated magnets the 
summability of the series is weak, leading to inaccurate or ambiguous resultaj'EHI. A second, and more fundamental, 
problem is the question of the matching between the predictions of the two field theoretical perturbative approaches: 
how can the existence of a non trivial fixed point around two dimensions be reconciled with the absence of such a 
fixed point around four dimensions ? Such a problem is, in fact, common to several field theories, relevant for type II 
superconductors, electroweak-phase transitions, smectic liquid cristal, He3, etc. In this sense, the situation observed in 
the vectorial 0(N) model where a non trivial fixed point exists everywhere between two and four dimensions appears 
to be the exception rather than the rule. Moreover, in the case of frustrated systems one encounters a more striking 
problem: there is an incompatibility between the symmetries of the NLcr and LGW models. Indeed, as said above, at 
the fixed point the symmetry of the NLcr model action is enlarged to 0(4), a phenomenon that cannot occur within 
the LGW perturbative approach. This raises serious doubts on the perturbative analysis of the LGW model away 
from d = 4. Reciprocally, the perturbative analysis of the NLcr model, based on a Goldstone mode expansion, predicts 
an 0(4)/0(3) fixed point, in contradiction with the perturbative LGW results and the experimental and numerical 
situation in d = 3. All this suggests that nonperturbative features could play a major role and thus imposes to go 
beyond the standard perturbative approaches. 



2 



III. EXACT RENORMALIZATION GROUP APPROACH 



To clarify this situation we have employed an approach based on the concept of effective average actionliillid, F/. [</>], 
which is a coarse grained free energy where only fluctuations with momenta q > k have been integrated out. The 
field cj> stands here for an average order parameter at the - running - scale k, the analog of a magnetization at this 
scale. In the limit k — > 0, all modes have been taken into account and I\. identifies with the usual effective action T. 
At the scale k = A - the overall cut-off - no fluctuation has been integrated-apd T k identifies with the microscopic 
action. The fc-dependence of I\. is controlled by an exact evolution equationliaiiJ: 



— --Tri(T k + R k ) — 



Tr{(T^+R k )-^\ , (4) 



where t = In k/ A. The trace has to be understood as a momentum integral as well as a summation over internal 
indices. In Eq.(|4|), 1^ is the exact field- dependent inverse propagator - i.e. the second derivative of T k - and R k is 
the effective infrared cut-off which suppresses the propagation of modes with momenta q < k. A convenient cut-off 
is provided by: Rk(q) = Z k q 2 / (exp(g 2 j k 2 ) — I), where Z k is the field renormalization. The effective average action 
is a functional invariant under the symmetry group of the system and thus includes all powers of all invariants - 
and their derivatives - built from the average order parameter. Thus, Eq.(^|) is a nonlinear functional equation, too 
difficult to be solved exactly in general. Therefore T k must be truncated. One - fruitful - possibility, is to perform an 
expansion of Y k around the minimum in order to keep a finite number of monomials, in the invariants while including 
the most relevant derivative termsl3. We have considered the following truncationta: 

r fe = J d d x | ^ v</> Q6 Wa& + y McaVM 2 + y ( ~ - « fc ) 2 + ^ A (5) 

where {cJfc, /z&, Zk) are the coupling constants that parametrize the evolution of the model with the scale, 

while p = Tr t <jxj> and r = ^Tr^ffxp) 2 — |(Tr t (j)(j)) 2 are the two independent O(N) (g> 0(2) invariants built from the 
average order parameter. Note that, since we have considered the generalization of the initial Heisenberg spin system 
to ./V-component spins, the order parameter consists in a real N x 2 matrix <\> a \,. Apart from the "current term" 
(eab^caV^cb) 2 - the terms in (^) are those appearing in the usual LGW action that realizes the symmetry breaking 
scheme of frustrated magnets. Indeed for X k and fik > 0, the minimum of the action is realized by a configuration of 
the form </>^ m — \/KkRab, where i? a & is a matrix built with two orthonormal iV-component vectors. The symmetry 
of this minimum is a product of a diagonal 0(2) group and a residual 0(N — 2) group. The symmetry breaking 
scheme is thus O(N) ® 0(2) -> 0(N - 2) ® 0(2) diag . Note that for <f> ab = 0™ fc m one has: p = 2n k and r = so 
that Eq.(|^) corresponds indeed to a quartic expansion around the minimum. At this minimum, the spectrum of the 
model, consists in 2N — 3 Goldstone modes and three massive modes: one singlet of mass mi = n k \ k and one doublet 
of mass 77i2 = Kfe/ifc which correspond to fluctuations of the relative angle and of the norms of the two vectors 4>i and 
(j)2- Let us finally stress that, without the current term, the truncation Eq.(||) is however not sufficient in our case. 
This term plays a crucial role since, for N = 3, it allows the model to enlarge its symmetry from 0(3) (£> 0(2) to 
0(3) (g> 0(3) ~ 0(4) at the fixed point, above d = 2. The current term is systematically discarded in the perturbative 
treatment of the LGW model around four dimensions, for the - correct - reason that it is power-counting irrelevant. 
Here we can include it in our ansatz since it is anyway present in the full effective action F^. In fact, we must include 
it since it becomes relevant somewhere between two and four dimensions, the formalism we use being in charge to 
decide where it is important. 

We have derived the flow equations for the different coupling constants n k , X k , Mfcj w fe an d Z k using Eq.(^) and 
Eq.(||)E3. The explicit recursion equations for the dirpensionless renormalized coupling constants g\ are too lengthy 
to be displayed here and not particularly illuminatingjij. Their general structure is given by: 

^ = Dgl+F i [{g r j },{T j [m^]}} (6) 

where D stand for the dimension - including the anomalous dimension - of the coupling g\. In Eq.(|^), the functions 
Fi involve the whole set of coupling constants {gj}- More important for our purpose is that the flow equations 
also involve the whole spectrum - Goldstone as well as massive excitations - of the theory. These masses occur as 
the arguments of some "threshold" functions {Tj} that govern the coupling or decoupling of the massive degrees 
of freedom. This is precisely this aspect that allows the method to go beyond the pure Goldstone mode expansion 
performed with the NLct model around d = 2. 

In practice, we have studied the fixed point structure of the flow equations (^) everywhere between two and four 
dimensions. Immediately around two dimensions, the fixed point masses 777* r and m^ 7 " diverge as 1/e. This means a 
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decoupling of the excitations associated to the fluctuations of the modulus and of the relative angle of <\>\ and cj>2- The 
physics is, as expected, fully controlled by the Goldstone modes. The flow equations (||) then degenerate to those of 
the O(N) ® 0(2) -> 0(N - 2) ® 0(2) d . iag NLcr model at one loop. For N = 3, they admit a fixed point for which the 
model is 0(4) symmetric, attesting the validity of the perturbative NLcr model approach. However, the flow equations 
actually admit another - unstable - fixed point, absent in the NLcr model approach, which moves toward the stable 
one while d is increased. At a critical dimension d c ~ 2.87, the two fixed points collapse and disappear. Above d c 
no other stable fixed point has been found, in agreement with the LGW approach around four dimensions. Note also 
that we have checked that, around four dimensions, an expansion of our equations in powers of the coupling constants 
allows to recover exactly the one-loop RG equations found within the LGW approach. The transition in d = 3 is 
thus, rigorously speaking, first order. However, the proximity of d c with d = 3 lets open the possibility of a very 
weak first order phase transition in this dimension. We have checked that there indeed exists, in d — 3, a range of 
parameters for which the RG flow is very slow. As a consequence the correlation length, although finite, is very large 
and a pseudo-scaling behaviour - an "almost second order phase transition "EJ - can be observed. Note however that 
£ remaining always finite, a true universality should not be expected. This could explain both the observation of a 
scaling behaviour as well as the existence of a broad spectrum of effective critical exponents in experiments and in 
numerical simulations. To test this scenario, we have computed effective exponents by linearizing the flow equations 
around the minimum, taken into account the </) 6 -like terms in our ansatz. We have checked that the inclusion of 
monomials of order six in the fields is sufficient to ensure the stability of our results. We have found: v = 0.53, 
7 = 1.03 and /3 — 0.28, which lie in between the various sets of exponents found experimentally and numerically. At 
this stage, one has to mention a work of Pelissetto et al. who have studied the O(N) <g> 0(2) — * 0(N — 2) §ljQ(2)diag 
LGW model at six-loop order in fixed dimension d = 3 and have reached a somewhat different conclusion!!!]. From 
their analysis, they have concluded that there exists a non trivial fixed poipt in the Heisenberg case. This is in contrast 
with a previous three-loop calculation in which no fixed point was foundBEa. The question is thus that of the stability 
of the perturbative results with respect to the number of loops or to the Borel resummation procedures. However, 
amazingly, the exponents found by Pelissetto et al. - v = 0.55(3), 7 = 1.06(5) and [3 = 0.30(2) - are very close to our 
pseudo-exponents, making the question of the existence of a true fixed point somewhat academic. Finally, we have 
found a true stable fixed point in d = 3 for N larger than a critical value N c (d = 3) ~ 5. For N = _£. using </> 10 -like 
terms, we get v = 0.707 and 7 = 1.377 that lie in between the error bars of the Monte Carlo datalid which display 
a clear second order behaviour with v = 0.700(11) and 7 = 1.383(36). Note that no result has been reported within 
the six-loop perturbative calculation in this case. 

Using an exact renormalization group approach, we have reached a global understanding of the physics of frustrated 
Heisenberg magnets including a matching between the previous perturbative predictions performed around the upper 
and lower critical dimensions and a good agreement with experimental and numerical data for the critical physics in 
d=3. 
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